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<N ■ Abstract 



We study, both experimentally and numerically, the Anderson localization phenomenon in tor- 
sional waves of a disordered elastic rod, which consists of a cylinder with randomly spaced notches. 
We find that the normal-mode wave amplitudes are exponentially localized as occurs in disordered 

Q ■ solids. The localization length is measured using these wave amplitudes and it is shown to decrease 

as a function of frequency. The normal-mode spectrum is also measured as well as computed, so its 
level statistics can be analyzed. Fitting the nearest-neighbor spacing distribution a level repulsion 
parameter is defined that also varies with frequency. The localization length can then be expressed 
as a function of the repulsion parameter. There exists a range in which the localization length 

,—1- \ is a linear function of the repulsion parameter, which is consistent with Random Matrix Theory. 

Tij- ■ However, at low values of the repulsion parameter the linear dependence does not hold. 

rn 
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FIG. 1: One-dimensional rod used to measure localization. In the experiment, the number of 
rods N = 50, R = 1.28 cm, d = (di) = 7.2 cm, and e = 1.016 mm; the values A = 0.35, 
rj = 0.65 were used. The largest frequency considered is less than 100 kHz so the lowest wavelength 
Amin = T^— > 2-R and the system behaves indeed as ID. The value c = 3140 — was measured for 

/max -5 

torsional waves in the aluminum alloy we have used. 

The Anderson localization phenomenon is a very important subject in condensed matter 
physics since it is crucial to understand the transport properties of materials. As a matter of 
fact, the original work of Anderson [1] is among the most cited papers in twentieth century 
physics and it is at the core of many papers, not only in solid state studies but also in optics. 



cold atomic gases, microwaves and acoustics 



12| . The theory of Anderson localization 



studies the alterations brought about on the localization of the electronic wave functions 
by disorder in the system. In a perfect lattice a band spectrum arises with extended wave 
functions for the allowed energy levels. However, if the system presents random imperfec- 
tions, for example the presence of strange atoms in an otherwise perfect structure or when 
there are unit cells of different size, wave functions can be localized, affecting the transport 
properties of the system. 

Anderson localization can also be studied using elastic vibrating systems, such as the 
quasi one-dimensional rod shown in Fig. [T] The system consists of A^ rods of radius R with 
lengths di, i = 1, . . . , N, joined by smaller cylinders of length e <t^ di,Wi, and radius r = rjR, 
where the coupling constant rj is such that < r] < 1. According to the nature of the family 
of numbers {di} different phenomena are observed. We should remark that we measure not 
only the normal-mode frequencies but also the wave amplitudes. In this sense, the analysis 
of elastic vibrations is more complete than what can be done in the quantum-mechanical or 
optical cases, since for such systems wave functions cannot, in general, be observed. 

In this Letter we shall study torsional vibrations in a disordered rod by taking the family 
{di} as a set of uncorrelated random numbers with a uniform distribution in the interval 



[d{l — A), d{l + A)], where d = {di) is the average of di and A measures the disorder. This 
allows us to study the Anderson localization phenomenon in elastic systems; in particular, 
we can observe how the wave amplitudes decay exponentially with a localization length ^. 
To perform the measurements we used the electromagnetic acoustic transducer (EMAT) 

n 

developed by us [13|. The EMAT consists of a permanent magnet and a coil, and can 
be used either to detect or excite the oscillations. The transducer operates through the 
interaction of eddy currents in the metallic rod with a permanent magnetic field. According 
to the relative position of the magnet and the coil, the EMAT can either excite or detect 
selectively compressional, torsional or flexural vibrations. Used as a detector, the EMAT 
measures acceleration. The experimental setup is described in detail in Ref. [141]. This 
transducer has the advantage of operating without mechanical contact with the rod. This 
is crucial to avoid perturbing the shape of the localized wave amplitudes. 

Before comparing the experimental data with the theoretical calculations, we shall con- 
sider what we will call an independent rod model [IJ]: The small rods of length di are 
independent from one another when 7] — )■ 0. In this case, the z-th rod is excited when the 
driving force has a frequency / equal to fn = ^, where c is the speed of torsional waves 
and n is an integer number. The other rods are in general not excited since dj for j ^ i is 
usually different from di. The amplitude of the vibration decreases when one moves away 
from the i-th rod so the wave is localized. It could also happen that some other length, say 
dj, can be almost equal to di. In this case the amplitude could then present two maxima. 
The first case is shown, for example, in Figs, ^bl) and|5|^c), while in Figs. MJo) and|5|^d) the 
second case is apparent. The independent rod model therefore provides a qualitative argu- 
ment to understand why all normal modes of a disordered rod are localized. Furthermore, 
when the disorder is very small, that is, when A ^ 1, the same argument shows that almost 
all the rods can be excited with a driving force of frequency / ~ ^ and the localization 
length ^ grows and could even exceed the total length of the complete rod; the amplitudes 
are then extended. 

If the cylinders have random lengths the elastic vibrations are waves on a random struc- 
ture, which is analogous to what happens with the Schrodinger wave functions in a random 
potential, frequency playing the role of energy. The independent rod model then shows that 
introducing disorder in {di} is a way to simulate diagonal disorder in a quantum mechanical 
one-dimensional tight-binding Hamiltonian, where the coupling r/ between nearest neighbors 



is a constant [15|. In this case the Anderson locahzation phenomenon also occurs. 

We shall now compute and measure the localization length ^ for a disordered rod. This 
quantity can be defined in at least three ways [16]: using the exponential decay of the trans- 
mission coefficient; by the Lyapunov exponent of the transfer matrix, or by the exponential 
decay of the wave amplitude envelope. Because the wave functions are in general not acces- 
sible experimentally, the first two definitions are normally used in the literature. However, 
we do have access to the wave functions, so we shall use the last definition. 

To obtain the experimental localization length ^ from the wave amplitudes, as those 
shown in Fig. [21 we consider a frequency interval and measure all the amplitudes with an 
eigenfrequency within this interval. An exponential is fitted by a least-squares procedure to 

the amplitude envelope 

(\x — X* I \ 
7 ) (1) 

at X > X* where ^4 is a constant. Here x* and ^ are the position of the maximum and the 

localization length of the wave amplitude, respectively. When the wave amplitude presents 

more than one maxima, only the highest one is used in the least-squares procedure. In this 

way the localization length was obtained as a function of the normal-mode frequency. To 

calculate the corresponding numerical wave amplitudes the Poincare map method [17'] was 

used and ^ was obtained in a similar fashion as before. It should be remarked that, as has 

been the case in all numerical calculations we have previously compared with experimental 

results, an effective value of the coupling parameter r] is necessary (See Ref. [l3| where this 

is justified). The calculations and the experiment coincide extremely well with each other 

using only this adjustable parameter. 

The localization length ^ as a function of frequency is given in Fig. [31 One should men- 
tion the overall agreement between theory (squares) and experiment (dots). The numerical 
values are obtained from an ensemble of 5000 families {rfj} but the experimental values were 
measured only for 50 eigenfunctions of a single rod. It is to be noted that ^ decreases with 
frequency. 

To obtain the wave amplitudes, such as those shown in Fig. [H the spectrum of the disor- 
dered rod must first be obtained. This is the case both numerically and in the laboratory. 
We are then provided with an extra bonus: the statistical properties of the elastic spectra 
which render themselves to studies like those analyzed in random matrix physics and in 



quantum chaos 18|, [l9|. In what follows we shall provide the nearest-neighbor spacing dis- 
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FIG. 2: Examples of experimental wave amplitudes for the disordered rod of Fig. [H The localization 
length ^ is extracted from the envelopes of these plots of the amplitudes as a function of the position. 

tribution p{si), where si = /f^^Zf\ is the normalized spacing, and show how the distribution 
varies as a function of frequency. 

We first consider an ensemble of 5000 disordered rods from a numerical point of view. 
The richness of the elastic ID system shows itself off when trying to fit the nearest-neighbor 
spacing distribution along the spectrum. We see that a more general function is needed than 
those calculated for the Gaussian orthogonal (GOE), the Gaussian unitary (GUE) and the 
Gaussian symplectic (GSE) random matrix ensembles. As a matter of fact, to obtain the 
level repulsion parameter a, the distribution 
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FIG. 3: (Color online) Average of the localization length (squares) as a function of the normal-mode 
frequency. The frequency average was done in windows of 20kHz. The experimental localization 
length for a single rod is given by the dots. It should be noted that ^ decreases as the frequency / 
grows. 

and A, B are constants such that the distribution and the average of the spacing are normal- 
ized (20|, was fitted to the numerical data by means of a least-squares procedure. Examples 
of numerical Pa{s) for different regions of the frequency spectrum are shown in Fig. |H The 
repulsion parameter a varies with the frequency. This is observed also in the experimental 
values of a shown in Fig. |5l 

Given these results one could also consider the relationship of ^ with a. This has been 
done when dealing with a disordered quasi-lD wire in mesoscopic physics. The Dorokhov- 
Mello-Pereyra-Kumar (DMPK) equation was derived [21| as well as a relationship between 
the localization length, the elastic mean free path l^i and the Dyson parameter /3 = 1,2,4, 
which characterize the Gaussian orthogonal, unitary and sympletic ensembles, respectively. 
The Dyson parameter (5 is also associated to the nearest-neighbor spacing distribution of 
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FIG. 4: (Color online) The nearest-neighbor spacing distribution p{s), obtained for an ensemble of 
5000 rods and average over a window frequency interval of 20 kHz, changes with frequency. From 
left to right the center of the interval is 28 kHz, 38 kHz and 78 kHz, respectively. 
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FIG. 5: (Color online) The nearest-neighbor spacing distribution p{s) obtained experimentally, 
with an average over a window frequency interval of 20 kHz, for a typical rod varies with the 
frequency. From left to right the center of the interval is 20 kHz, 24 kHz and 48 kHz, respectively. 



the quantum energy spectrum. Indeed a = P for 13 = 1,2, 4. It was shown 22] that a linear 
relationship between ^ and /3 exists 



e=[(iV,-l)/3 + 2]4 



(4) 



where Nc is the total number of channels (for a review see 23J). Eq. (j4]) was verified for 
/3 = 1,2,4 in magnetoconductance experiments 24l |. 

The plot C, vs a for the elastic rods is given in Fig. [61 The points were obtained for the 
ensemble of rods. For a = 1,2,4 this repulsion parameter coincides with /?; one can see 
from the crosses in this figure that indeed ^ is a linear function of a. But this is not true 
at low values of a, where a higher order polynomial fits the theoretical values better. The 
experimental results for our rod also show this, which is a new effect, not obtained with the 
usual random matrix ensembles. 
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FIG. 6: (Color online) Average localization length ^ (squares) as a function of the repulsion 
parameter a. The values a = 1, 2, 4 (crosses) correspond to the random matrix theory predictions, 
with Nc = 2, ^e = 0.115 m, for which ^ is a linear function of a (dashed line). For a < 1 this 
dependence is no longer linear. The dots correspond to the experimental values obtained for a 
single rod. 

A qualitative explanation for the decrease of ^ when a — )■ is again provided by the 
independent rod model. When 77 — > 0, the frequency spectrum, defined by the independent 
random sequence {dj}, is also random. As is well known in random matrix theory the nearest- 
neighbor spacing resembles then a Poisson distribution, that is, the repulsion parameter 
tends to zero. On the other hand, ^ also diminishes in this case, since the wave amplitudes 
are localized in a single small rod. We therefore understand, from a qualitative point of 
view, the structure of the curve shown in Fig. |6l 

To conclude, in this Letter we measured the exponential localization of elastic wave 
amplitudes in disordered rods as well as the spectrum and calculated the level repulsion 
parameter a. It was found that the localization length presents two regimes. For a > 1 
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the localization length grows linearly as a function of the level repulsion parameter. This 
agrees with the predictions of the DMPK equation for the Dyson random matrix ensembles. 
However, there is a region, a < 1, in which the localization length is not a linear function 
of a. 

This work is supported by DGAPA-UNAM under project IN111311. RAMS was sup- 
ported by CONACYT under project 79613. 
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